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We  have  devised  a fully  empirical  prtKcdurc  for  generating  chemically  accurate  (/.<■.,  within  a few 
kcal  mol-‘)  potential  hypersurfaces  for  triatomic  bimolecular  reaction  systems.  The  procedure  is 
described  and  applied  to  the  H,  system.  Calibration  of  the  H3  empirical  surface  by  comparison  with 
the  best  ah  initio  H 3 surface  suggests  considerable  utility  of  the  empirical  procedure. 

1,  INTRODUCTION 

Refined  and  accurate  ah  initio  potential  energy  hypersurfaces  can  now  be  obtained 
for  simple  {i.e.,  few  nuclei  and  few  electrons)  chemically  reactive  systems  such  as 
Hj,  FH2,  and  HLi2.‘  Certainly,  the  ah  initio  route  is  preferable  to  all  other  computa- 
tional schemes  in  cases  demanding  modest  expense  {c.g.,  less  than  $10  000  in  1976 
U.S.  dollars)  for  the  several  hundred  hypersurface  points  required  to  yield  a good  fit  to 
chemically  interesting  regions  of  the  surface.  Nevertheless,  the  popularity  of  semi- 
empirical  surfaces  of  the  LEPS  (London-Eyring-Polanyi-Sato)  and  DIM  (Diatomics- 
In-Molecules)  varieties  attests  to  the  continuing  utility  of  inexpensive  (albeit  of  inde- 
terminate accuracy)  surfaces  in  practical  contexts  (e.g.,  quasiclassical  trajectory  studies 
of  reaction  dynamics). 

We  present  here  a tully  empirical  scheme  for  the  computation  of  potential  hyper- 
surfaces for  triatomic  bimolecular  reaction  systems.  The  following  sections  describe 
our  empirical  approach,  calibrate  its  accuracy  by  comparison  to  the  best  available 
H3  ah  initio  ground  potential  hypersurface,  and  comment  upon  extensions  and 
limitations  of  the  method. 

2.  RESUME  OF  THE  MODEL 

Empiricism  precludes  theoretical  rigour.  In  the  empirical  limit,  intuition  replaces 
formal  and  numerical  certainty.  Our  empirical  model  builds  upon  the  set  of  intuitive 
assumptions  used  in  the  Johnston-Parr  BEBO  (Bond  Energy-Bond  Order)  method 
for  the  construction  of  collinear  reaction  coordinates  of  H-atom  transfer  systems: 

Assumption  1 — The  amount  of  “ chemical  bonding  ” is  constant  at  each  point  along 
the  reaction  coordinate. 

Assumption  2 — The  total  interaction  potential  is  the  separable  sum  of  three  two-body 
potentials. 
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Assumption  3 — Empirical  correlations  (Pauling’s  rule,  Badger’s  rule,  etc.)  can  be  used 
to  calculate  bond  energies  and  other  molecular  properties  from  structural  and  spectro* 
scopic  data. 

A.  BEBO  COLLINEAR  REACTION  COORDINATE**^ 

In  symbols,  for  the  reaction: 

A+BC  - A-^“  B C AB+C 

where  B is  a hydrogen  atom  (although  modified  BEBO  treatments  have  been  developed 
for  non  H-atom  transfers)***  and  /Ixy  is  the  “ bond  order  ” of  diatom  XY,  we  have: 

'Iab+Wbc^I.  (1) 

independent  of  position  along  the  reaction  coordinate  (wxy  = i for  each  diatom  well 
separated  from  its  corresponding  atom) 

1^TOTAl(^AB>  ^Bc)  = 1^Ab(^Ab)+ l^BcC^Bc)"!' ^Ac(^Ac)»  (2) 

with  the  energy  referenced  to  the  A+BC  limit  and  with  /?ac  = -^ab+^bc  for  the 
collinear  system 

/?xv(^xy)  ““  -^e.xY  — “0.26  In  Wxy  and  (^^) 

^Xy('*Xy)  = ■f^e.XY^f*'^  (=^Xy) 

where  and  Z),  are  the  isolated  diatom  equilibrium  internuclear  separation  and 
dissociation  energy,  R(n)  and  D(n)  are  the  dressed  diatom  structure  and  dissociation 
energy  for  bond  order  n,  and  p is  a “ bond  index  ” (close  to  unity)  giving  the  func- 
tional dependence  of  bond  energy  upon  bond  order. 

Using  valence  bond  concepts,  Johnston  and  Parr*  evaluated  (2)  by  taking  K^b  and 
l^Bc  to  be  attractive  potentials  [calculated  from  eqn  (36)]  and  to  be  a repulsive 
potential  [represented  by  a scaled  anti-Morse  fit  to  the  Hirschfelder-Linnett  calcula- 
tion* of  the  H2*(6*ZJ)  interaction  energy].  They*  chose  to  calculate  the  bond  index 
from 

_ 0.26  In  jPJe)  ... 

^ Rm-R. 

where  e and  are  the  potential  parameters  (well-depth  and  position  of  the  minimum) 
for  corresponding  rare  gas  diatoms  (e.g.,  HCi  is  the  rare  gas  pair  which  is  used  to  mimic 
a low  bond  order  version  of  dressed  H2).  Thus,  the  BEBO  potential  energy  along  the 
collinear  reaction  coordinate  is: 

V{n)  = Z)..2(l-/i‘’*)-Z)..i(l-n)'’'  + Z)..3B(n-n*)“**^[l+^(/i-«*)®*‘*']  (2') 

where  n = /Ibo 

B = 0.5  exp  [~fiAR^], 

P is  the  Morse  potential  parameter  [=(fc*/2D,)‘'*  where  A:,  is  the  force  constant]  for 
diatom  AC, 

AR,  = 

and  subscripts  1,  2,  3,  are  used  in  place  of  subscripts  AB,  BC,  AC. 

By  comparison  with  ab  initio  results,  Truhlar*  showed  that  the  original  BEBO  model 
(using  p = 1.041  for  H2  and  p = 1.036  for  HF)*  generated  amazingly  accurate  reac- 
tion coordinates  and  saddle  point  heights  for  H3  and  FH2  systems.  We  update  his 
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H3  conclusions  slightly  in  fig.  1 and  2.  Ab  initio  results*  given  in  each  figure  are  the 
most  accurate  collinear  Hj  computations  presently  known,  with  stated  uncertainty 
limits  of  lg;|  kcal  mol"‘  in  surface  point  energies  and  f5xl0"“  A in  the 
coordinate  location  of  the  saddle  point.  The  BEBO  reaction  coordinate  (/.e.,  the 
minimum  energy  path),  defined  by  eqn  (1)  and  (3o)  is  even  more  accurate  than  previ- 
ously recognised''  by  comparison  with  an  earlier  //n7/o  computation;’  the  maximum 
deviation  (at  the  saddle  point)  is  now  ~().01  A {cf.  fig.  1).  This  result  is  not  subject  to 


^,/A 


Fig.  1 . — BEBO  ( ) reaction  coordinate  for  ground-state  H3.  Ab  initio  points  ( x ) are  from  ref.  (8). 


Fig.  2.— BEBO  potential  energy  profiles  for  ground-state  Hj.  Parameters  for  old  BEBO  ( ) and 

new  BEBO  I ( ) and  II  ( ) plots  are  described  in  Appendix  B.  Ab  initio  ( x ) points  are  from 

ref.  (8). 

updating  [except  [jerhaps  through  a revision  of  the  constant  appearing  in  Pauling's 
rule,  eqn  (3a)-]  The  old  BEBO  form  of  the  potential  energy  along  the  reaction  co- 
ordinate [from  eqn  (2'),  using  = 1.041]  is  also  strikingly  accurate  when  compared 
with  Bowen  Liu’s  ah  initio  potential  energy  profile*  {cf.  fig.  2),  but  this  agreement  must 
be  viewed  with  suspicion.  Fig.  2 shows  new  BEBO  potential  energy  profiles  which 
incorporate  new  potential  parameters  for  Hcj  derived  from  scattering  experiments'® 
leading  to  = 0.997  together  with  the  ab  initio  H2*(6^Sj;)  repulsive  potential." 


130  POTENTIAL  HYPERSURFACE  FOR  BIMOLECULAR  REACTION  SYSTEMS 

The  updated  parameterization  of  p and  V^^c,  designated  “ new  BEBO  I ” in  fig.  2, 
gives  too  low  a saddle  point  height.  However,  slight  alteration  of  for  Hei  from 
the  best  scattering  potential  value  of  2.963  A to  an  adjusted  value  of  2.93  A increases 
the  bond  index  p to  a value  of  1.012  and  raises  the  collinear  saddle  point  to  close 
agreement  with  the  ah  initio  result  (cf.  curve  designated  “ new  BEBO  II  ” in  fig.  2). 
These  (sensitive)  changes  illustrate;  (a)  that  the  success  of  the  old  BEBO  H3  collinear 
potential  energy  profile  (based  upon  crude  Lennard-Jones  potential  constants  for 
He^  plus  an  inaccurate  H2  repulsive  potential)  is  partly  fortuitous,  and  (b)  that  accurate 
potential  data  still  yield  quite  reasonable  saddle  point  heights  and  potential  energy 
profiles. 

We  conclude  that  the  original  Johnston-Parr  BEBO  modeF  and  its  updated  ver- 
sions provide  a fine  representation  for  the  H3  reaction  coordinate.  Numerous  other 
systems  have  been  treated  by  BEBO  methods,^"*'’’*^"*^  with  the  result  that  calculated 
saddle  point  heights  are  often  within  1-2  Real  mol"^  of  experimentally  observed 
Arrhenius  activation  energies.'* 

B.  EXTENSION  OF  BEBO  TO  FULL  POTENTIAL  HYPERSURFACES 

We  retain  the  BEBO  picture  of  a triatomic  reaction  system  A B C as  a 

pair  of  “ dressed  " diatomic  molecules  AB  and  BC  with  repulsion  between  end  atoms 
A and  C.  At  any  point  along  the  reaction  coordinate  and  in  the  vicinity  of  atom  C, 
diatom  AB  has  a bond  energy  that  is  diminished  relative  to  the  isolated  molecule;  the 
corresponding  dressing  of  diatom  BC  occurs  due  to  the  presence  of  atom  A.  Eqn 
(3fi)  and  (2')  simply  evaluate  the  altered  diatom  bond  dissociation  energies.  For  H3 
(and  most  other  reaction  examples),  a bond  index  close  to  unity  signifies  that  the  total 
attractive  interaction  remains  nearly  constant;  hence,  the  repulsive  interaction  is  the 
dominant  contribution  to  the  increase  in  total  energy  along  the  reaction  coordinate. 
However,  if  we  leave  the  reaction  coordinate,  there  are  other  increases  in  total  energy 
that  we  view  as  simply  the  increase  in  potential  energy  due  to  the  displacement  of 
dressed  diatoms  AB  and  BC  from  their  equilibrium  internuclear  separations.  For 
example,  fig.  3 shows  the  BEBO  H3  reaction  coordinate  and  a point  Ri,  Ri  displaced 
from  the  minimum  energy  path.  Two  reference  points  (7?, a,  Ri^  and  R^,,  R^^,)  along 
the  reaction  coordinate  are  also  indicated.  With  respect  to  point  ^1.,  Ri^,  the  dis- 


Fio.  3. — Reference  point  notation  for  extended  BEBO  calculations.  The  solid  line  segment 
is  the  BEBO  reaction  coordinate  for  ground-state  Hj. 
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placed  point  Ri,  refers  to  the  same  dressed  diatom  AB  that  was  used  to  construct 
the  BEBO  reaction  coordinate  point  /?ja,  /?2a  (since  the  AB  internuclear  separation  is 
unchanged)  but  a stretched  dressed  diatom  BC.  With  respect  to  point  R^,,  Rzb,  the 
displaced  point  /?,,  Rz  refers  to  two  stretched  diatoms  (AB  and  BC)  since  both  AB  and 
BC  internuclear  separations  have  altered.  There  is  an  infinity  of  reference  points 
{/?!,,  Rzi)  along  the  reaction  coordinate  and  the  total  energy  depends  upon  the  degree 
of  stretching  and/or  compression  associated  with  each  choice.  For  all  our  calcula- 
tions, we  choose  the  intuitively  satisfying  reference  point,  viz.,  the  one  that  minimizes 
the  total  energy:*" 

d y(Ru  ^ ^ 

d/7i 

where/;,  = 1 — defines  the  bond  order  and  hence  the  location  of  the  reference  point. 

In  the  following  development,  we  neglect  the  repulsive  interaction  (which  is 
not  a function  of  the  bond  order,  since  it  depends  solely  on  /J^);  for  any  geometry: 


B— C 


we  have: 

/?,  ^ R^^  'rlR^Rz  cos  (6) 

and  is  determined  either  from  a scaled  anti-Morse  function’-*  or  from  our  fit 

to  the  ab  initio  H2*(/>^It)  potential.'*  The  total  energy  excluding  repulsion  is: 

F(/?„  R,)  Ml  ! Mz  (7) 

where  M,  /)?{! —exp[—/iy(/ii  — /??)]}-  (8) 

/y?  = ikU2D<iy''‘  (9) 

and  A-y.,  — exp(y,  f (10) 


Reference  diatom  quantities  D",  R",  and  /;,  are  evaluated  from  eqn  (3)  and  (4).  The 
Herschbach-Laurie  version*’  of  Badger’s  rule,'*  eqn  (10),  is  used  to  calculate  dressed 
reference  diatom  force  constants  A?_,;  the  parameters  5’,  and  rl>|  are  determined  from  a 
fit  to  experimental  data  on  an  inter-related  set  of  molecules  or  several  electronic  states 
of  the  same  molecule.'’ 

The  iterative  solution  of  eqn  (5)  is  discussed  in  Appendix  A.  Once  the  reference 
value  «i  is  determined  from  eqn  (5),  F(R,,  Rz)  is  calculated  from  eqn  (3),  (4)  and  (7). 
The  total  energy  is  the  sum  f"Ac(/?3)  1 y(R„  Rz)- 


3.  EXAMPLE:  H,  GROUND  POTENTIAL  HYPERSURFACE 

Fig.  4 shows  equipotentials  for  Hj  at  = O”.  30°,  60°  and  90°  calculated  by  our 
empirical  procedure  (as  described  in  Appendix  B).  As  0 increases,  the  saddle  point 
height  increases  and  the  minimum  energy  path  shifts  out  to  larger  internuclear 
separations;  hence  the  collinear  reaction  pathway  is  preferred.  Within  the  context 
of  BEBO  empiricism,  collinear  surfaces  have  the  lowest  saddle  point  energies  due  to 
the  purely  repulsive  character  of  the  AC  interaction.  Close  inspection  of  the  collinear 
surface  in  fig.  4 reveals  that  the  coordinate  location  of  the  saddle  point  is  slightly 
shifted  outward  (to  7?,  = 7?^  = 0.956  A)  relative  to  the  Pauling  rule’"  result  [from 
eqn  (3a),  /?,  = = 0.922  A].  This  shifting  is  due  to  the  relative  steepness  of 

attractive  (F^b,  Fbc)  and  repulsive  (F,^c)  interactions.  For  H3,  the  anti-Morse 
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0 50  100"  150  .Too  \ 250  0 50  l”00  150  , 200  2 50 

^,/A  R:  /A 

Fig.  4. — Empirical  potential  hypersurfaces  (0  = O’,  30’,  60’  and  90  ’)  for  ground-state  Hj. 


repulsion  listed  in  Appendix  B rises  steeply  as  /?j  decreases,  while  the  Morse  attrac- 
tions arc  quite  modest  since  c;  MJi’Ri  ~ 0 in  the  vicinity  of 

It  is  important  to  obtain  an  assessment  of  the  quantitative  success  of  the  empirical 


:o  zo  30 

#'i3  J 


> Fio.  5. — Potential  energy  profiles  along  the  minimum  energy  paths  of  empirical  H3  ground  hyper- 
f surfaces. 

surface  we  have  constructed ; this  assessment  can  best  be  accomplished  by  comparison 
of  both  its  static  and  dynamic  characteristics  relative  to  those  of  the  best  H3  ah  initio 
surface.®  One  static  calibration  involves  the  shape  of  minimum  energy  paths.  Fig.  5 
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shows  potential  energy  profiles  along  the  minimum  energy  paths  for  Hj  0 - 0”,  30"’, 
60°  and  90°  empirical  surfaces;  these  profiles  compare  favourably  with  accurate 
ah  initio  results.^*  The  bending  potential  is  faulty,  however  (e.g.,  the  saddle  point 
height  fort?  = 90°  is  only  23  kcal  mol"'  compared  with  approximately  30  kcal  mol"' 
for  an  accurate  ah  initio  surface  “). 

A multitude  of  empirical,  semi-empirical,  and  ah  initio  H , surfaces  [rev  iewed  in  ref. 
(23)]  have  been  calculated.  We  here  restrict  our  discussion  to  comparison  with  a few 
widely-used  surfaces  (LEPS,  DIM,’'*  Porter-Karplus,^*  etc.).  The  calibration  is  al- 
ways Liu’s  ah  initio  collincar*  and  non-ci>llinear“  surfaces,  as  defined  by  an  array  of 
137  collinear  points  and  surface  fit"  together  with  1 15  non-collinear  points “ and  cor- 
responding surface  fits.^'’  We  choose  to  compare  the  points  themselves;  table  1 
shows  standard  deviations  of  the  energies  of  various  surface  points  as  a function  of 


Table.  1— Standard  deviations  of  empirical  and  semi-empirical  sori ace  points  irom 

LIU  CALIBRATION  POINTS. 


standard  deviations  for; 

surface 

reference 

F - 20/ 

y-  40/ 

kcal  mol"' 

kcal  mol  ' 

kcal  mol"' 

Yates-Lester* 

21 

1.15 

1.28 

1.48 

Kafri-Berry 

this  work 

1.52 

2.20 

2.82 

improved  LEPS 

h 

1.12 

1.17 

1.38 

Porter-Karplus, 

surface  2'^ 

25 

1.8.2 

2.38 

2.63 

‘ This  surface  is  intended  to  be  an 

accurate  fit  to  Liu 

calibration  points. 

^ London -Eyring- 

Polanyi  Sato  form  incorporating  iih  initio  repulsion"  plus  modified  Morse  attraction  (see  Appendix 
B);  Sato  parameter  adjusted  to  0.025  8 to  yield  the  correct  saddle  point  height.  *■  The  most  popular 
version  of  this  potential.” 

maximum  potential  energy  sampled.  For  L 20,  30  and  40  kcal  mol"',  samples  of 
145,  175  and  214  ah  initio  points"-^^  arc  available.  Probably,  only  points  below  a 
maximum  potential  energy  of  20  kcal  mol"  * arc  important  for  thermal  reaction  pro- 
ces.ses.  We  are  carrying  out  quasiclassical  trajectory  computations  to  gauge  the 
dynamical  accuracy  of  various  H3  surfaces  at  this  time,  we  note  that  our  empirical 
surface  is  simple  and  inexpensive  to  use  in  such  dynamical  computations. 

In  contrast  to  LEPS  and  Porter-Karplus  surfaces  which  contain  adju.stable  para- 
meters, the  present  empirical  surface  (based  on  old  BEBO^-^)  is  not  adjustable. 
However,  new  BEBO  I and  II  reaction  coordinates  (see  fig.  2 and  Appendix  B)  which 
incorporate  much  better  defined  potential  parameters  apparently  do  not  yield 
better  surfaces;  for  example,  standard  deviations  are  about  50"o  higher  for  ’’  ri  20, 
30  and  40  kcal  mol"'  calibration  points  compared  with  old  BEBO  results. 

Another  popular  surface  which,  in  principle,  is  not  adjustable  is  DIM  (Dia- 
tomics-In-Molecules).^'*  Recent  extended  DIM  results^"  show  that  the  method  is 
stable  to  basis  set  variation  and  is  reasonably  accurate  provided  that  an  ad  hoc 
neglect  of  overlap  factors  is  adopted.  The  best  agreement  of  extended  DIM  was 
found  for  cases  1(a)  and  11(a)  [see  ref.  (28)]  which  predicted  collinear  saddle  point 
heights  of  12.5  and  12.0  kcal  mol"'.  We  can  anticipate  that  standard  deviations  of 
DIM  points  from  Liu  ah  initio  points*'^  will  be  somewhat  higher  than  the  surfaces 
considered  in  table  I. 

Standard  deviations  may  be  misleading,  since  they  are  sensitive  to  both 
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coordinate  and  energy  errors  (e.g.,  a trial  surface  isomorphous  to  the  true  surface, 
but  shifted  significantly  in  coordinate  space,  would  display  large  standard  deviations; 
reaction  dynamics  computations  using  such  an  isomorphous  surface  would, 
however,  probably  be  quite  accurate).  In  lieu  of  dynamical  computations,  a better 
measure  of  surface  suitability  for  treatment  of  reactive  processes  might  be  the  accur- 
acy of  the  isotropic  (i.c.,  spherically  symmetric)  part  of  the  potential,  now  known  ex- 
perimentally for  Fig.  6 shows  the  experimental  Fq  together  with  Porter- 

Karplus,  Yates-Lester  and  our  ^'o  values;^®  we  seem  to  be  in  quite  close  accord  with 
experiment. 


4.  SOME  PROJECTIONS  AND  CONCLUSIONS 

We  have  empirically  treated  various  other  systems  (FH2,  CNH2,  etc.)  with  as  much 
apparent  success  as  for  Hj.  The  old  BEBO^'^  prescription  for  the  collinear  reaction  co- 
ordinate seems  to  guarantee  reasonable  saddle  point  heights  and  geometries, 
even  though  the  treatment  of  the  triplet  repulsion  between  end  atoms  is  dubious.  We 
liave  also  been  able  to  extend  our  empirical  surface  to  4-atom  systems  such  as  H2  + D2, 
but  calibration  of  the  results  is  difficult  for  these  more  complex  cases. 


Air 


#/A 


Fki.  6. — Spherically  symmetric  (i.e.,  isotropic)  part  of  the  H: 

identified  in  table  I ; experimental  data  arc  from  ref.  (29).  


Yates-Lester, Porter-Karplus. 


H potential.  .Sources  of  data  are 
experimental, Kafri-Berry, 


Extensions  of  the  empirical  approach  described  in  this  paper  to  non-reactive  and 
excited  potential  hypersurfaces  appear  feasible.  At  present,  we  cannot  recommend  a 
good  procedure  for  dealing  with  multi-valent  atoms  (e.g.,  oxygen  atom  in  water)  and/or 
systems  with  preferred  non-collinear  paths. 

The  greatest  utility  of  our  empiricism  may  well  be  for  the  treatment  of  reactive 
collisions  of  simple  triatomic  systems,  but  this  application  remains  to  be  tested  and 
documented. 

We  wish  to  thank  Drs.  Bowen  Liu,  George  Schatz,  Don  Fleller  and  B.  H.  Choi  and 
Profs.  Phil  Certain,  K.  T.  Tang  and  Don  Truhlar  for  discussions  and  materials. 
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APPENDIX  A: 

SOLUTION  OF  EQUATION  (5) 

In  order  to  solve  eqn  (5)  to  find  and,  hence,  the  reference  point  along  the  reaction 
coordinate  (via  eqn  (3a)],  we  take  partial  derivatives: 


dK(/?.,  R^) 
dn. 


. /?2)_^r-gP?  dMjdR^i  dM,8D°,  SM, 
1 /=!  _ c*rii  drii  dfi° 


(IK 

'c;R°i  c*w,  ‘ cD^  tn,  )_ 


From  eqn  (3),  (7),  (8),  (9)  and  (10),  we  have: 

dOydrii  = piD°Jni  (A2) 

dXfjen,  = -2^fD?exp[-/?f(/?i-/??)]{l-exp[-/??(/?i-/??)]}  (A3) 

dRljdni  = —0.26lrii  (A4) 

dMJdDi  = MJD°1  (A5) 

dMjep°,  = 2{R,-R^)Dlexp[~p^(R-Rl)]{\~exp[-p1{R^-R°,)]}  (A6) 

BPV8kl,=p°,(2kl,r^  (A7) 

dklJdR°,  = d,kl,  (A8) 

a/??/aD?  = -p\{2DX)-\  (A9). 

Corresponding  partial  derivatives  with  respect  to  «2  are  obtained  via  eqn  (1)  which 
yields  the  ratio:  driPdrii  = — 1. 

The  explicit  form  of  eqn  (Al)  is: 

AJriy-Ailn^  = 0 (AlO) 

where  /4,  = /j,(A/,— D?)+T,X,  (All) 

r,  -0.26(2D?y??)  + D?A?(/?,-/??)(0.26^,-/j,)  (A12) 

and  X,  =exp[-^?(/?, -/??)](! -exp[-/??(7?, -7??)]}.  (A13) 

Newton-Raphson  iteration^’  was  used  to  solve  eqn  (Al): 

"j  + i = nj—F(8FI8rii)-^  (A14) 

where 

F=8V(RMI8ni  (A  15) 

8FI8n,  = (Q,-A,)lnl+{Qz-A2)lnl  (A16) 

and  Qi  = p,At+T,K, (0.26St-pi)l2+p,TtKt+0.26p’lK, {0.263, -p,)D<i 

-2DM0.26+(Ri-m0.263,-p,)/2]}^cxp[-p-,{R,-R<t)]{l-2cxp[~p^XRi-m)- 

(A17). 

To  start  the  iteration,  we  use  the  initial  condition: 

nu  „.r,  = n?/(n?+nS)  (A18) 

with  n°,  = exp[(/?,.,— /?,)/0.26].  (A19) 

Convergence  (n^+i— ^ 10"*)  is  typically  reached  after  3 iterations. 
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APPENDIX  B: 

II3  SURFACE  PARAMETERS 

Three  sets  of  H3  surfaces  were  calculated  by  using  old  BEBO,  new  BEBO  I,  and 
new  BEBO  11  collinear  reaction  coordinates  (c/.  fig.  2).  In  all  sets,  we  took:  £>?  = 
109.5  kcal  mol-‘,  0.741  2 A,  and  A:?  , = 286 exp(4.00-0.25/??)  kcal  mol"* 

l or  old  BEBO,  p,  ^ 1.041  and  \\c  ^ o’.25Z)5{exp[-2y95(/?3-/?3)]  + 2exp[-/?5(/?3- 
/?3)]1-  Tor  new  BEBO  1.  Pi  - 0.996  and  V xc  is  the  ab  initio  repulsion."  For  new 
BEBO  11,  p^  — 1.012  with  the  same  ah  initio  repulsion.  Fig.  4 shows  old  BEBO 
computations. 


For  most  surface  calculations,  we  employed  Morse  functions  for  My  and  Mi  in 
eqn  (7).  However,  the  Morse  function  for  HiiX^Zg)  is  a poor  fit  to  the  ab  initio 
potential^’  at  R < and  even  deviates  significantly  (by  more  than  4 kcal  mol“‘) 
at  long  range  ( ~2  A).  For  the  thermal  H + H2  reaction  (i.e.,  for  surface  points  below 
20  kcal  mol"'),  the  inadequacies  of  the  H2  Morse  function  are  not  serious,  but  higher 
energy  processes  require  a better  potential  fit. 

We  have  devised  a modified  Morse  function  for  Hi{X^l.'^)  which  is  accurate  to 
better  than  0.5  kcal  mol"'  at  R > 0.57  A: 

My  = Z)'J{l-exp[-/^"(/?,-/?^)(l+Z.,)]}^ 

Ly  = h^c\rp[—hy{Ry—R°y-hiY] 

/>o  = 0.019  959  95 
by  - -2.217  032  9 
^2  = 0.512  704  12 
b^  = 0.283  162  54 
by  = 0.836  707  93 
hi  = 2.251  996  4. 

All  three  H3  surfaces  (old  BEBO,  new  BEBO  1,  and  new  BEBO  II)  use  the  modified 
Morse  potential.  Solution  of  eqn  (Al)  requires  a slight  alteration  of  working  equa- 
tions given  in  Appendix  A since  My  has  a different  form.  We  shall  present  details  of 
the  modified  Morse  fitting  procedure  and  hypersurfacc  calculations  elsewhere. 


where 
At  R R°: 

while  at  R~^  R'\'. 


(Bl) 

(B2). 
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